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New theoretical calculations are performed to investigate the Coulomb proximity and angular momentum ef- 
fects on multifragmentation picture for **Kr+!!*:'4Sn collisions at an incident beam energy of 35 MeV/nucleon. 
Charge and isotopic distributions and the mean neutron-to-proton ratios of the fragments are reproduced within 
the microcanonical Markov chain calculations on the basis of Statistical Multifragmentation Model. It is shown 
that the Coulomb interactions and angular momentum effects are very important to reproduce isotopic com- 
position of nuclear fragments in peripheral heavy-ion collisions at Fermi energies. Our results imply that it is 
possible to investigate in laboratories the modification of structure parameters of fragments, such as the symme- 
try energy coefficient, at subnuclear densities in dense environment of other species. 
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I. INTRODUCTION 


During the peripheral heavy-ion collisions at Fermi 
energies (20-50 MeV per nucleon), a considerable amount 
of angular momentum could be transferred from the interac- 
tion region to the excited projectile and target residual nu- 
clei, and this can lead to significant changes in their multi- 
fragmentation [1—6]. Additional long-range forces caused by 
the complicated Coulomb interaction between the target and 
projectile-like sources are essentially involved in the process 
[3, 7]. The multifragmentation in the presence of the external 
Coulomb field offers a possibility to study, experimentally, 
the effects of this long-range force, which are very important 
for disintegration of matter [2]. This is also necessary for 
construction of a reliable equation of state (EoS) of nuclear 
matter at subnuclear densities. Another motivation of these 
studies is that similar conditions for nuclear matter take place 
during the collapse and explosion of massive stars and in the 
crust of neutron stars [8, 9], where the Coulomb interactions 
of dense electron environment change the fragmentation pic- 
ture. It is generally assumed that the statistical equilibrium 
regarding the fragment composition at subnuclear densities 
should be established in these astrophysical cases. There- 
fore, the analysis of the observables obtained in laboratory 
experiments with statistical models is a proper way to obtain 
knowledge on stellar matter. Previous studies of isospin com- 
position of the produced fragments were found to be espe- 
cially important for determining the strength of the symme- 
try energy during fragment formation in hot and diluted en- 
vironments [10-15]. In central heavy-ion collisions at Fermi 
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energies (20-50 MeV per nucleon), high excitation and high 
densities can be reached [16]. They become a suitable tool 
to study EoS of hot nuclear matter and the nuclear liquid-gas 
phase transitions at subnuclear densities. As shown previ- 
ously, one can study the properties of hot fragments in the 
vicinity of other nuclear species by means of multifragmenta- 
tion [15]. The angular momentum effect is usually neglected 
in this case, since the impact parameters are small. 


As demonstrated by several studies [12, 17—21], using sta- 
tistical ensemble approach within Statistical Multifragmen- 
tation Model (SMM) [22], charge and isotope yields, frag- 
ment multiplicities and temperatures, and correlations of var- 
ious fragment properties were successfully analyzed using 
ALADIN data. This was also achieved in the analysis of 
the experimental data [23] obtained at the MSU laboratory 
at 50 MeV/nucleon [24—26], and in the analysis of TAMU 
data [27, 28]. In these studies, the symmetry energy of frag- 
ments was one of the main model parameters governing the 
mean (N)/Z values, the isoscaling parameters, and the iso- 
topic composition of the fragments. In our theoretical inter- 
pretations [12, 24-26] for ALADIN and MSU experiments, 
we considered the formation and decay of single thermalized 
source and the averaged Coulomb interaction of fragments 
(the Wigner-Seitz approximation), since a direct positioning 
of fragments in the freeze-out volume has minor influence on 
their charge and isotope distributions. We have justified the 
assumption of the formation and decay of single thermalized 
sources for relativistic peripheral collisions, and central colli- 
sions of heavy nuclei around the Fermi energy with our suc- 
cessful predictions for ALADIN and MSU data. We believe 
that it is possible to extract important information on multi- 
fragmentation and properties of fragments in peripheral colli- 
sions at Fermi energies as well. The new fragment partitions 
can be obtained by including the Coulomb effects caused by 
the proximity of colliding target and projectile nuclei, as well 
as the effects by large angular momentum transfer to the mul- 
tifragmenting sources. We believe that a long-range Coulomb 
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interaction of the target- and projectile-like sources changes 
the fragmentation pattern and leads to a predominant midra- 
pidity (neck—like) emission of light and intermediate mass 
fragments (IMF, with charge numbers Z =3-20). In Refs. 
[28, 29], such experiments have already been analyzed with 
statistical models. However, there were no systematic the- 
oretical investigations of the Coulomb and angular momen- 
tum effects on multifragmentation picture at these reactions, 
especially on the isotope yields which are crucial for astro- 
physical applications. The angular momentum may lead to 
more neutron-rich IMF production and to anisotropic emis- 
sion with respect to the projectile and target sources, as it was 
reported in Refs. [3, 7]. In this paper, we theoretically in- 
vestigate the influence of angular momentum and Coulomb 
interactions on the charge yields, the neutron to proton ra- 
tios of particles for peripheral **Kr+!!7!74Sn collisions at 
35 MeV per nucleon. This is a quite typical reaction, and our 
selection is partly motivated by recent FAZIA experiments 
[30]. Our investigation on interpretation of FAZIA experi- 
ment is underway. For the simulation of the reactions, we 
consider the break-up of a single source **Kr in the prox- 
imity of a secondary source !!7:!?4Sn. The calculations are 
carried out within the Markov chain version of the statis- 
tical multifragmentation model (SMM), which is designed 
for a microcanonical simulation of the decay modes of nu- 
clear sources [3, 31]. This method is based on producing the 
Markov chain of partitions which characterize the whole sta- 
tistical ensemble. In this method the individual fragment par- 
titions and coordinate positions of fragments in the freeze-out 
volume are generated. They are selected by the Metropolis 
algorithm and we can take into account the influences of an- 
gular momentum and Coulomb interactions for each spatial 
configuration of primary fragments in the freeze-out volume, 
similar to Refs. [1, 2]. 


Il. STATISTICAL METHOD FOR 
MULTIFRAGMENTATION CALCULATIONS 


In the microcanonical SMM, a statistical equilibrium at low 
density freeze-out region is a basic assumption. Nucleons and 
nuclear fragments are included in the breakup channels, and 
the laws of conservation of energy Ex, momentum, angular 
momentum, mass number A and charge number Z are consid- 
ered. The breakup channels and the compound-nucleus chan- 
nels are also included, and competition between all channels 
is permitted. The SMM includes the conventional evaporation 
and fission processes occurring at low excitation energy as 
well as the transition region between the low and high energy 
de-excitation regimes. In the thermodynamic limit, SMM 
is consistent with liquid-gas phase transitions when the liq- 
uid phase is represented by infinite nuclear clusters [32], that 
allow connections for the astrophysical cases [6, 33]. The sta- 
tistical weights of all breakup channels partitioning the sys- 
tem into various species are calculated. The decay channels 
are generated by Monte Carlo method according to their sta- 
tistical weights. In the Markov chain SMM [3, 31], we use 
ingredients taken from the standard SMM version developed 
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in [22, 34, 35] which was successfully used for comparison 
with various experimental data: Light fragments with mass 
number A < 4 and charge number Z < 2 are considered as el- 
ementary particles with the corresponding spins (nuclear gas) 
that have translational degrees of freedom. The fragments 
with mass number A > 4 are treated as heated nuclear liquid 
drops. Therefore it is possible to study the nuclear liquid-gas 
coexistence in the freeze-out volume. Free energies F4 z of 
each fragment are parameterized as a sum of the bulk, surface, 
Coulomb and symmetry energy contributions 
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The bulk contribution is given by FẸ; = (Wo — T? /£o)A, 
where T is the temperature, the parameter £ọ is related to the 
level density, and Wọ = 16 MeV is the binding energy of in- 
finite nuclear matter. Contribution of the surface energy is 
l; = ByA*?((T2-T7)/(T2+T7)P4, where By = 18 MeV is 
the surface energy term, and Te = 18 MeV is the critical tem- 
perature of the infinite nuclear matter. In the standard SMM 
version the Coulomb energy contribution is E$ , = cZ*/A'”’, 
where c denotes the Coulomb parameter obtained in the 
Wigner-Seitz approximation, c = (3/5)(e7/ro)(1 — (p/po)!3), 
with the charge unit e, rọ = 1.17 fm, and pọ is the nor- 
mal nuclear matter density (0.15 fm~°). However, within this 
Markov-chain SMM we directly calculate the Coulomb in- 
teraction of non-overlapping fragments in the freeze-out by 
taking into account their real coordinate positions. The sym- 
metry term is Ee = y(A — 2Z)/A, where y = 25 MeV is 
the symmetry energy parameter. All the parameters given 
above are taken from the Bethe-Weizsaecker formula and cor- 
respond to the assumption of isolated fragments with nor- 
mal density unless their modifications in the hot and dense 
freeze-out configuration follow the analysis of experimental 
data. Since our previous analysis [13] confirms the trend of 
a decreasing symmetry energy as one approaches conditions 
comparable to the multifragmentation regime in agreement 
with previous findings [11, 12, 15, 24, 26, 28], we have used 
y = 14MeV for Z < 9, y = 16MeV for Z = 10 — 17, and 
y = 18-19 MeV for Z > 18. In this work, we use p = po /6 to 
freeze out density for better evaluation of Coulomb and angu- 
lar momentum effects. Usually, we generate about 10° Monte 
Carlo events to provide sufficient statistics. 


HI. CALCULATIONS WITH ANGULAR MOMENTUM 
AND COULOMB INTERACTIONS 


In this study, we investigate peripheral nucleus-nucleus 
collisions at 35MeV/nucleon with corresponding relative 
velocities of the projectile and target around 45-70 mm/ns. 
At the initial dynamical stage of such a collision, the pro- 
jectile nucleons interact with target nucleons and some en- 
ergetic products of this interaction can leave the nuclei as 
pre-equilibrium particles. The kinetic energy of colliding nu- 
clei can also be converted into the excitation energy of pro- 
jectile and target residues. Therefore, the relative velocity 
between the residues decreases as well. These excited target 
and projectile-like sources decay afterwards. 
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Table 1. Values for the sources assumed to be formed according to 
the peripheral collisions for **Kr+!'*:!"4Sn reactions. 


E,(AMeV) Z, <A,('?Sn) A, (Sn) Weight ang. mom. 
2 34 77 82 0.13 30 
3 33 75 80 0.19 30 
4 32 7 71 0.32 40 
5 30 68 73 0.25 50 
6 29 66 70 0.08 60 
7 28 64 68 0.03 70 


The projectile and target-like sources will not be far from 
each other before disintegration since nuclear multifragmen- 
tation is a fast process within a characteristic time around 
100 fm/c. We believe that at these short distances the long 
range Coulomb field of one of the sources affects the break 
up of the other one. In this way, we can deal with the mul- 
tifragmentation in a double nuclear system, since it is a new 
physical situation with respect to the standard multifragmen- 
tation of a single isolated source. 

According to our estimates from the energy conservation, 
their relative velocity should decrease to ~60 mm/ns, at an 
excitation energy around 5 MeV/nucleon transferred to the 
residues. In this case, the first and second source will be sep- 
arated by ~15fm in a time of 100 fm/c. The decay of the 
two excited sources in such a double system is determined 
by the short-range nuclear forces. However, the presence of 
an external Coulomb field (for each source) may influence 
the composition of the produced fragments and their relative 
positions. Especially, an additional Coulomb barrier will pre- 
vent disintegration of the sources into many small pieces. It 
should be noted that during the evolution of a double system 
we must take into account its total center of mass conserva- 
tion without a separate constraint in the freeze-out volumes of 
disintegrating sources. On the other hand, we include the an- 
gular momenta (rotation) of the separate sources, which can 
be transferred after the collision. It will also influence the 
positions and sizes of the fragments at the freeze-out [1, 3, 7]. 

We show the results for multifragmentation of the 
projectile-like source (the first source **Kr) by assuming the 
Coulomb field coming from the center of the target source 
(the second source !!7:!74Sn). We have assumed that the first 
source flies along the Y-axis, and the second one flies in the 
opposite direction (related to the center of mass of the dou- 
ble system). This separation axis may slightly deviate from 
the initial beam axis. The location of the second source is 
taken as Ry= -15fm and Rz=15 fm with respect to the first 
source. The peripheral collision is assumed to take place in 
the Y — Z plane. The coordinates in Z-axis are determined 
by the sizes of colliding nuclei, as well as by their possi- 
ble repulsion after the collision. The angular momentum 
axis is assumed to be X-axis. We believe that this relative 
space configuration of the sources is suitable for investigat- 
ing the Coulomb and angular momentum effects. The pre- 
equilibrium emission of few nucleons during the dynamical 
stage may decrease the excitation energy and relative veloc- 
ity of the residues. One may take into account these effects 
in the statistical approach by changing the corresponding in- 
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put and by using the ensemble of the sources [12, 13] with 
adequate parameters. We also know from many theoretical 
and experimental works [17, 22, 36], that the relative yields 
of IMF do not depend much on the size of the sources in the 
multifragmentation regime. But in our previous work [37], 
we also investigated the behavior of sources of the same size 
and isospin content as the colliding nuclei (**Kr). 
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Fig. 1. The charge yield of fragments, in the cases of without 


(full circles) and with angular momentum (open circles), after mul- 
tifragmentation of the projectile “Kr source. These sources are 
assumed to be formed in the peripheral **Kr + !'7:!*4Sn collision at 
35 MeV/nucleon, and its disintegration is affected by the Coulomb 
field of the target source as shown in Table 1. 


A. Charge and isotope distributions, and neutron to proton 
ratios 


We have investigated the angular momentum and Coulomb 
field influences on the charge and isospin contents of the 
produced fragments. It is important to analyze the new char- 
acteristics of fragment distributions, which are crucial for in- 
terpretation of many experiments on heavy-ion collisions at 
Fermi energies. After the break-up of the sources we cal- 
culate the Coulomb propagation of produced hot fragments 
by taking into account the Coulomb interactions of particles 
in the double system. In order to clarify the modification of 
the multifragmentation picture caused by the new effects and 


S$20507-3 


Nihal Buyukcizmeci et al. 


1.4 oles | a a agi at ins ] Lae a ae ] a ae a |i sel a | ae fe | J ETI ] ae ian a | rey 
1.3 É 84k r+ 4Sn K 
i —@— with angular momentum | 
: F —O— without angular momentum 7+ 
N L J 
à lee =| 
Vv i 4 
1c am 
4k r+ "Sn 
Z=3-9 y=14 MeV 
N Z=10-17 y=16 MeV 
N 
< Z=18-20 y=19 MeV 


Fig. 2. The neutron-to-proton ratio (N)/Z of fragments produced at 
the freeze-out density p = po/6. 
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Fig. 3. Predicted isotope distributions for Z = 6,12 and 18 frag- 
ments for **Kr+!!?Sn (top panel) %*Kr+!*Sn (bottom panel) colli- 
sions. The solid lines correspond to the present results with angular 
momentum, the dashed lines show results without angular momen- 
tum. For our predictions y = 14, 16 and 19 MeV values are used 
for Z=6, 12 and 18 isotopes, respectively. These y values are in 
agreement with our previous predictions in [12, 13, 26, 37]. 


compare it with experimental data in the future, we apply the 
secondary de-excitation of the hot fragments which can lead 
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to important consequences especially for the isospin compo- 
sition of final fragments. 


For our calculations we made following assumptions: We 
take the sources obtained from the **Kr projectile at different 
excitation energies, 2, 3, 4, 5, 6 and 7 MeV/n, as shown in 
Table 1. Liquid-gas phase transition theory is applied in this 
energy interval to explain nuclear multifragmentation on the 
basis of SMM. At FE, < 2 MeV/nucleon the compound nu- 
cleus and fission channels are found to be dominant. We have 
taken two sources with the same charge and the same excita- 
tion energy, but with two different (NV)/Z ratios corresponding 
to those of **Kr+!!*Sn and 84Kr+!4Sn, as in experiment [30]. 
Then, firstly, we made SMM calculations for each source and 
found the charge and isotope distributions and (N}/Z ratios of 
fragments in the presence of second source (!!*!74Sn) without 
and with angular momentum. Angular momentum values are 
randomly selected with an increasing order (see Table 1). Af- 
terwards, we have taken mixture of all sources with weights, 
as shown in Table 1, corresponding to their excitation ener- 
gies, which are related to the impact parameters. After this 
mixture is determined, we propose that we can obtain similar 
conditions to compare with experiments. 


In Fig. 1, we show the charge yields of cold fragments with 
and without angular momentum conservation. Angular mo- 
mentum values are selected as shown in Table 1. As seen in 
Fig. 1, an angular momentum favors emission of large nearly 
symmetric fragments (like a nuclear fission) since the system 
in the freeze-out needs to have a large moment of inertia in 
order to minimize the rotational energy and to maximize the 
entropy. It is in a competition with the second source through 
the Coulomb interaction which prevents to emit an IMF with 
a large charge number. 


The initial value of neutron-to-proton ratio of the projectile 
source **Kr is 1.33 while 1.28 and 1.42 are corresponding to 
initial (N)/Z ratios of °*Kr+!!?Sn and 84Kr+!4Sn. In Fig. 2, 
we see that angular momentum leads to increasing (N)/Z val- 
ues of light IMFs in the case of strongly asymmetric decay. 
By the increasing moment of inertia of the system, it favors 
a bigger phase space of the reaction [3]. This trend may be 
responsible for many isospin observables. 


We show in Fig. 3 the variation of isotopic distribution for 
Z=6, 12 and 18 for 8*Kr+!*4Sn and **Kr+!!?Sn reactions. We 
have demonstrated that isotopic distributions are very sensi- 
tive to Coulomb and angular momentum effects for peripheral 
reactions at Fermi energies. On the basis of our findings, we 
belive that it should be taken into account these new effects 
for the realistic description of experimental data, e.g., FAZIA 
[30]. Possible modifications for the symmetry energy term 
can also be investigated by means of the isotopic distributions 
of projectile fragments. In previous works, it is seen that the 
calculations for Z = 6 (Z < 10 small fragments) agree well 
with the data at y = 14 (as in MSU and ALADIN analysis 
[24, 26]), for Z = 12 at y = 16 while the calculations for 
Z = 18 with y ~ 19 [13]. In a forthcoming paper, we will 
introduce the results of detailed analyses and interpretation of 
FAZIA [30] experimental data, on the basis of present calcu- 
lations. 
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To verify our new-found trends we also performed the 
similar calculations for lighter and heavier systems, e.g., 
84Kr+*4Kr and !?’Au+!?’Au collisions. In all cases, we 
have got the similar qualitative modifications of the standard 
multifragmentation picture related to the effects of angular 
momentum and the Coulomb field of sources. 


IV. CONCLUSION 


As a result, we have theoretically investigated the charge 
and isotope distributions and (N}/Z values of fragments after 
the multifragmention of the Kr-like projectiles in peripheral 
84K r+!!2,124Sn collisions around the Fermi energy within the 
microcanonical Markov chain approach on the basis of the 
statistical multifragmentation model. Coulomb and angular 
momentum effects originated after the collision dynamics are 
taken into account for the first time in this study. We demon- 
strated that conservation of angular momentum and compli- 
cated Coulomb interactions caused by the proximity of target 
and projectile-like sources in the freeze-out stage produce sig- 
nificant changes in the multifragmentation picture. There ap- 
pears to be new fragment formation trends, such as an asym- 
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